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Abstract 

In this paper, the effects of a newly presented structure of 
unified power flow controller (UPFC) which is based on two 
shunt voltage-source converters and a series capacitor, are 
investigated on the dynamic stability of power system. A 
power system including the newly presented structure of 
UPFC is introduced to examine this. The non-linear dynamic 
model of the system under investigation and its non-linear 
equations are presented as well. The linear equations are 
achieved by linearizing the non-linear equations around the 
operation point. A modified Heffron-Philips model is 
proposed through the linearized equations of the 
investigated system. The state equations of the system are 
then extracted using the modified Heffron-Philips model. 
The damping controller of UPFC is then presented and its 
control parameters are calculated by phase compensation 
method. Finally, the most effective input signal of UPFC 
control signals is selected using some techniques. The 
validity of the theoretical investigations is confirmed by 
simulation results obtained in MATLAB/Simulink. 
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Nomenclature 


M 

Inertia constant (second) 

Xe, M g = 1,2 

Transformer g reactance (pu) 

X sek for k = a,b,c 

Equivalent reactance of capacitor in k phase 

(pu) 

X. for i = d,q 

Steady state synchronous reactance of the 
generator on i axis (pu) 

K 

Direct axis transient synchronous reactance of 
the generator (pu) 

x bv 

Reactance of transmission line (pu) 

x „ 

Reactance of transformer (pu) 

Q. 

Dc link capacitor (pu) 

c„ 

Series capacitor (pu) 

R E s for g = 1,2 

Equivalent resistance of shunt converter g (pu) 

I b for i = d,q 

Current in infinite bus on i axis (pu) 


II 

J*, 

Current of phase k in transmission line on i 
axis (pu) 

I sei for i = d,q 

Current of phase k in series capacitor on i axis 
(PU) 

h, s) f or S = U 

i = d,q 

Current in converter g on axis (pu) 

u 

Voltage at dc link (pu) 

V se k for k = a, b,c 

Series capacitor voltage of k phase (pu) 

K 

Voltage at infinite bus (pu) 

v t : 

Voltage at generator terminal (pu) 

II 

£ 

Terminal voltage on i axis (pu) 


Reference voltage (pu) 

A T upfc 

Provided electrical torque by UPFC damping 
controller (pu) 

P 

Input mechanical power to the generator (pu) 

P 

Electrical power of the generator (pu) 

A 

State matrix 

B 

Input matrix 

C 

Output matrix 

D 

Direct transfer matrix 

X 

State vector 

Y 

Output vector 

U 

Input vector 

m 

Forward path of the UPFC damping controller 

K c (s),K„(s), 

K l (s),K m (s) 

UPFC damping controller transformer controls 

H(s) 

Transfer function of UPFC damping controller 

f 

Switching frequency (Hz) 

8 

Phase angle of System load (degree) 

A 8 

Phase angle deviation of System load (degree) 

S E , g for g = 1,2 

Phase angle of shunt converter g (degree) 

Adi,» for g = 1,2 

Phase angle deviation of shunt voltage source 
converter g (degree) 

m E, s for g = 1.2 

Pulse width modulation index of shunt 

converter g 

A m Ejs for g = 1,2 

Pulse width modulation index deviation of 
shunt voltage source converter g 

CO Q 

Fundamental angle frequency ( rad s 1 ) 

CO, 

Oscillation angle frequency ( rad s~' ) 

8(0 

Error signal ( rad V 1 ) 

D m 

Friction factor of generator 
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^UPFC 

Friction factor of the UPFC damping controller 

U k 

Input signal of the UPFC 

'V(M) 

System operation domain 

U, 

UPFC controller gain 

K a 

Gain of AVR in generator 

K , ~K n 

Constants of Heffron-Phillips model 

c,uk ’ p,uk ’ <7,u* ’ 

K^ K se,k 

Constants control of input signals 


Constant control of input signal 

T L 

Direct axis open circuit time constant of the 
generator (second) 

T a 

Time constant of the generator of AVR 

’Tu,T w 

Time constants of phase compensator 


Introduction 

Nowadays, considering the increasing growth of 
power systems in different countries and their mutual 
connections and also the great interest in transmitting 
the possible maximum power through the transmission 
lines, the maintenance and enhancement of power 
systems stability becomes more important [Shahir 
(2013), Kundor (1994)]. The modem power systems are 
designed to achieve high generation capacities, 
bringing more risk into power systems stability. One 
of the most important factors threatening power 
systems stability is the electromechanical oscillation of 
the generator. Dynamic instability is raised due to 
balance lose between the output electric power and 
input mechanical power, lack of damping torque, and 
sudden load distortion outbreak. The frequency of 
such oscillations falls in the range of one to few Hertz 
[Shahir (2013)]. Several factors such as operation 
condition, load characteristics, lines impedance and 
power, series capacitors, and voltage regulators 
influence this type of oscillations. 

One of the most effective solutions is to install power 
system stabilizer (PSS) in generator location to damp 
the electromechanical oscillations [Yassami (2010)]. 
The PSS generates a complementary signal to rapidly 
damp the distortions caused by oscillations of power 
system. The PSS usually considers frequency error 
signal as its input and generates a signal added to the 
voltage. The practical and simulation results show that 
the PSS will be highly effective in system oscillations 
damping if it is properly regulated. 

Recent technology development in semiconductor 
based devices and their control techniques field and 
high interest of power systems designers in utilizing 
flexible ac transmission systems (FACTS) for maximum 
power transmission has provided a new stability 


research field. The application of Applying FACTS 
devices brings about more reliability and less costs in 
different network arrangements under distortions 
occurrence conditions [Song (2004), Xia (2010)]. Among 
the FACTS devices, the voltage-source based converters 
considerably affect the transient state and improve 
dynamic stability of power systems [Hingorani (1999)]. 
The UPFC, one of the most important and effective 
FACTS device applied in series-shunt mode consists of 
a series and a shunt three-phase two-level converter 
connected to each other through a dc capacitor. The 
series converter controls the active and reactive powers 
of transmission line by controlling the amplitude and 
phase angle of series voltage and the shunt converter 
compensates the UPFC losses and stabilizes dc link 
voltage in addition to transmission line reactive power 
control [Sarvana (2009)]. Many different structure have 
been proposed for UPFC [Shahir (2012), Sadigh (2010)]. 
In [Sadigh (2010)], a new UPFC structure consisting of 
two shunt voltage-source converters and a dc capacitor 
installed between the mentioned converters is 
presented in which the UPFC parameters are controlled 
by the series capacitor. The voltage and the current of 
the series capacitor are controlled by the operation of 
the two shunt converters [Sadigh (2010)]. The control 
and the measuring systems of conventional UPFCs are 
separately designed due to the different insulation and 
protection processes of series and shunt converters. This 
leads to high time and cost expenditure of such devices 
design and manufacturing. In the structure of UPFC 
presented in [Sadigh (2010)], the evaluation bulk is 
decreased due to similar shunt converters application 
and as a result; the control and measuring systems 
design and manufacturing costs are significantly 
decreased. Also, in comparison with congenital UPFC, 
total VA rating for high active and reactive power flow 
control is decreased. In the structure of [Sadigh (2010)], 
the signals could be applied as UPFC input signals for 
stability and electromechanical oscillations control that 
is an advantage of that structure. Three signals can be 
applied for stability control in conventional UPFCs 
[Shahir (2012)], while four signals are used to control 
the operation of the newly presented UPFC. Reduction 
of the THD value of the injected series voltage and 
applying series capacitor as constant series capacitive 
compensator during converter 2 standby interval are 
other advantages of the structure presented in [Sadigh 
( 2010 )]. 

The modeling process of conventional UPFCs and its 
effect on dynamic stability is investigated in [Fang 
(2003)] and transient stability rate improvement is 
evaluated in [Sreenivasachar (2000)]. In this paper, the 
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effects of the UPFC structure presented in [Sadigh 
(2010)] on the dynamic stability and electromechanical 
oscillations damping of power system are investigated. 
The analysis starts by nonlinear equations extraction. 
The linear equations and the modified Heffron-Philips 
model are obtained by linearizing the nonlinear equations 
around the system operation point. The phase 
compensator technique [Larsen (1995), Chandrakar 
(2010)] is then applied and the proper controller is 
designed and regulated. Finally, the most effective 
input signal for UPFC control is obtained by simulating 
the system in MATLAB/Simulink to damp the 
oscillations of synchronous generator's rotor. 


Investigated System 


A sample power system including a single machine 
with infinite bus along with the UPFC presented in 
[Sadigh (2010)] is used to investigate the dynamic 
stability. The single line structure of the investigated 
power system is illustrated in Fig. 1. The utilized 
UPFC consists of two parallel transformers, a series 
capacitor and two shunt voltage-source converters 
controlled by pulse width modulation (PWM) technique 
[Shahir (2012)]. Here, R El and R E2 are neglected in 

comparison with X E l and X E 2 , respectively. 



FIG. 1 THE SINGLE LINE DIAGRAM OF INVESTIGATED SYSTEM 


Dynamic Modeling of the Investigated Power 
System 


Nonlinear Dynamic Model 


All reactance and resistances of the generator, 
transformers, and transmission line should be 
considered for nonlinear dynamic modeling of power 
system. Therefore, the following nonlinear equations 
are obtained assessing the system shown in Fig. 1: 


8 = a> 0 (a) - 1) 

P„, ~P„ - D, „ A w 


co = - 


M 


( 1 ) 

( 2 ) 


P, ~ E ‘, +E fJ 
^ q T > 

1 Art 


1 do 


Pfd ~ 


- Efd+ K a (V ref -V t ) 


V dc = 


3 m 


E,l 


4C, 


dc 


cos S E i sin S El 


Ie,u 


L E,\,q 


3 m 


4 C 


7 ^“ [cos 8 E2 sin S E2 ] 

'dc 


1 E,2,d 


1 E,2,q 


Ke = — f se,i f° r i = d,q 
^ se 

The followings are valid for (l)-(6): 

P e = V ud I t4 + V tq I tq -,E q =E' q+ (X d - X' d )I t4 ; 
V t =V t4 +jV Uq ,V ud =X q l tq ,V tq =E ' q -X' d I t4 ; 

d t,d = d se,d + 1 E\,d ’ If,q = d se,q + d E,\,q ' 


b \ m EoVdc . 
2 


I se,d = a x V b cos 8 + ^ " cos S E2 + c x E' q + 


, , d l m E,\ V dc 
2 


b im E ,2 V dc ■ x , d \ m E /dc • 


he, q = a \ V b sin 8 + ^-sin 8 E2 + 


sin 5 E ! 


r w X b 2 m E,2^dc x r , d 2 m E,\V d c 

I Eid = a 2 V b cosh + cosdp i +c- ) L„ + 


J E, 2 


hxa = a 2 V b sinh + 


b 2 m E,2 V dc . , d 2 m E,\Vdc . 

-smh r -,+ 


sin S E j 


, TS ■ X b 3 m E,2^dc „ r , d 2 m E,\^dc 

J E, 2 ,d = a } V b sill 8 + cos 8 E 2 + c 3 E q + 

T TS ■ X b 2 m E,2^dc ■ O d 3 m E,l^dc ■ x 

i E2q =a 2 r b smd-\ : sinh^j -1 1 sinh, 


E, 1 


I t,d = («1 + a 2 Wb cos 8 + {\+ b 2 ) E ^ ‘ /c cos S E 2 
+ (q + c 2 )E' q + (d x + d 2 ) E ± dc cos S E1 


(3) 

(4) 

(5) 

( 6 ) 


cos S E x 


cos j 


COS §j? j 


1 t, q = («i + a 2 )V b sin 8 + {b x + b 2 ) n ' E - 2 J dc s in S E 2 

,,1,1 , m E,iP dc . x 

+ {d x +d 2 ) — - — sin S El 

x bb = X E2 +X bv ;X al = X se +X El X le + X te X se ; 

X a 2 = X BB X E y,X a2 = X E 2 X E l ; X bl = X El + X te ; 

X b \ = X E \ + X te ;X b2 = X E2 X bv ;X b2 = X BB X se ; 

^ 1 =V, 2 + x bl x b2 ; X 2 = X iv + ; 

a e, i 

The other constants are presented in Appendix A. 
Linear Dynamic Model 

The linear dynamic model of the sample system is 
obtained linearizing (l)-(6) around the system operation 
point. The following equations are obtained after 
linearization: 
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AS = co 0 Aco 


A co = ■ 


A P„ - A P„ - DAco 


ae ;=- 


A£ /t/ = 


M 

-AE„+AE 


fd 


T' 

1 do 


-AE fd+ K a (AV ref -AV t ) 


(7) 

( 8 ) 

(9) 

( 10 ) 


AV dc = K 7 AS + K,AE' q - K 9 AV dc + K c EA 


Am 


E, 1 


+ ^c,SE,\^E,\ + ^-c,E,2^ m E,2 + ^c,SE,2^E,2 


( 11 ) 


AV se -K l0 AS + K xl E +K sl Am E 2 + K s2 Am E , 


+ K n AV dc + K s3 AS Ea + K s4 A8 EX 
The followings are valid in (7)-(12): 

A P e = K x A 8 + K 2 AE' q + K pEX Am EX + K pM - 


( 12 ) 


[A S E i 

E i (13) 


+ ^p,E,2^ m E,2 + ^p,SE,2^E,2 + ^ pd^dc 

A E q = K 4 AS + K 3 AE' q + K qE1 Am E ^ +K qMl A8 E3 

+ ^q,E,2^ m E,2 + ^q,SE,2^E,2 + ^qd^dc 

AV t =K 5 AS+K 6 AE’ q + K vEX Am E2 + K vSE4 AS E4 

+ K V 'E,2^ m E,2 + ^v,SE,2^E,2 + ^vd^dc 


(14) 


(15) 


The constants used in (11)-(15) are presented in details 
in Appendix B. 


The Proposed Modified Heffron-Philips Model 

Here, it is possible to propose the modified Heffron- 
Philips model of the system shown in Fig. 1 based on 
(7)-(12) to analyse and investigate the system stability 
against oscillations and to evaluate its small-signal 
stability [Shahir (2012), Kundor (1994)]. The proposed 
modified Heffron-Philips model is a linear model 
applied in low frequency oscillations and small signal 
stability investigations. Fig. 2 shows the block diagram 
of the proposed modified Heffron-Philips model 
resulted from (7)-(12). In this model, one of the 
linearized control signals of the sample power system 
including UPFC, A m E X , A m E 2 , A S El , and A S E2 is 
applied to the UPFC and controls the investigated 
UPFC's operation. 

The block diagram of proposed modified Heffron- 
Philips model shown in Fig. 2 contains 35 constants 
that can be defined as row vectors of K p uk , K q uk , 


K v ,uk > K CM k / and 

K k as follows to simplify the 

calculations: 





Kp,uk = \-Kp,E, 1 

K p,SE,l 

Kp,E, 2 

Kp,8E, 2] 

(16) 

Kq, u k = [K q ,E, 1 

Kq,SE,\ 

Kq,E,2 

KqJE, 2] 

(17) 

k vM = \.K vEX 

K v ,se, i 

K v ,E, 2 

K v ,SE, 2] 

(18) 

^ c.uk ~ \^c,E, 1 

K c ,SE, 1 

K c ,E, 2 

K-c,SE, 2] 

(19) 


K se,uk ~ 1 K s ,SE, 1 ^s,E, 2 ^s,SE,2 ] ( 20 ) 

The row vectors of (16)-(20), are applied on proposed 
modified Heffron-Philips system as the following 
input vector: 

U = [Am EX AcJ^j Xm E2 AS E2 \ (21) 



FIG. 2 THE PROPOSED MODIFIED HEFFRON-PHILIPS MODEL 
OF INVESTIGATED SYSTEM 


The State Equations of the Investigated System 

Based on (7)-(15), the dynamic model of the investigated 
system in the form of the state space equation can be 
expressed as follows: 

X = AX + BU (22) 

where X and U are defined as follows: 


X = [A S A co A E’ q A E fd AV dc AV se ] 


U = [Am 


E, 1 


A S El Am 


E, 2 


A5 E2 ] 


(23) 

(24) 


A and B matrices of (22) are considered as follows: 

(T 

o 


A = 


0 

co 0 

0 

0 

0 

h. 

0 

_ffi_ 

0 

K pd 

M 

M 

M 

k 4 

0 

k 3 

1 

K qd 

t;, 0 

T'dO 

T' 

1 d 0 

Tdo 

K a K 5 

0 

K a K 6 

1 

K a K vd 

T a 

T a 

T 

a 

T a 

Kq 

0 

K, 

0 

-k 9 

K 10 

0 

K n 

0 

K n 


(25) 
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0 

0 

0 

0 

Kp,E, 1 

Kp,SE, i 

Kp,E, 2 

Kp,SE, 2 

M 

M 

M 

M 

Kq,E,\ 

Kq,SE,l 

Kq,E, 2 

Kq.SE, 2 

T> 

1 d0 

T' 
1 d 0 

T< 

L d 0 

T' 

J d 0 

K a K v E> i 

K a K v gE' i 

K a K v ,E, 2 

K a K v ,sE,2 

T a 

T 

a 

T a 

T 

a 

K c,E,l 

K c ,se. i 

K c ,E,2 

K c ,SE, 2 

K s ,e, i 

K s ,se, i 

K s ,E,2 

K s ,SE, 2 


It is possible to design an appropriate controller using 
(22) and select the most proper input signal of UPFC 
from four input signals according to criteria detailed 
in the next section in order to maintain the system 
stability and to control the damping of power system 
against the rotor oscillations of investigated system's 
synchronous generator. Therefore, it is important to 
select and design appropriate input signal and 
desirable damping controller. 


UPFC Damping Controller 

The linearized model of the UPFC installed in power 
system explained in (22) can be shown in Fig. 3 [Wang 
(1997), Shahir (2012)]. The UPFC damping controller 
consists of two sections of input signal and UPFC 
controller. The input signal is applied on power 
system through the input transfer functions of K 0 (s ) , 
K l (s), and K c (s) . The state space equation obtained 
in (22) can be considered as follows [Wang (1997), 
Shahir (2012)]: 

0 o 0 0 AS 0 

-k m -d m -a d22 A co + -b d2 U (27) 

a d 31 a d32 fl d33_|L Z mJ V~^d2 _ 

"Ad" 

Y = \ c d\ c d2 c d2,\ A® +DU (28) 



FIG.3 THE CLOSED LOOP BLOCK DIAGRAM MODEL FOR UPFC 
DAMPING CONTROLLER 

In (27), k m and d m are the synchronizer and damper 



components, respectively, and z m consists of other 
state variables of (22). The oscillation angular 
frequency is obtained by ®, = yj co 0 k m . Here, the 
transfer functions presented in [Wang (1997), Shahir 
(2012)] are explained using (27) and (28). 

K c (s ) is the transfer function from U input to Ad) 
and obtained as follows: 

K c ( 5 ) = a d22 (si - a d33 ) b d 3 + b d2 (29) 

K 0 (s) is the function of transferring rotor speed to the 
y measuring point shown in block diagram of Fig. 3 
and denoted as follows: 

Ko( s ) = c d)( s I ~ a d33) 1 ~ + a d32 j + + c d2 j (^0) 

Kj ( s ) is the function of transferring U input to the y 
measuring point shown in Fig. 3 and calculated as 
follows: 

K rX s ) = C/3(‘ s/ -^ 33 )" V 3 + D (31) 

K m (s) is the transfer function of Ad to Ado and 
calculated as follows: 


K m ( s )-k m + dm +a d2i( s I a d3i'> \ a d2\ + a d)2 ( 32 ) 

®o l ®0 ) 

H(s), as a UPFC controller transformer, controls the 
operation of the system shown in Fig. 1 for different 
input signal values. In this investigation, the UPFC 
damping controller is considered ideal. An ideal 
UPFC damping controller should provide a pure 
positive electric torque in oscillation loop, damp 
electromechanical oscillations, and stabilize the 
generator rotor oscillations. Considering Fig. 3, the 
provided electric torque by UPFC damping controller 
could be calculated as follows [Wang (1997), Shahir 
(2012)]: 

(33) 

UPFC 1 ~K l (s)H(s) 

The torque magnitude provided by UPFC damping 
controller is expressed as follows: 

A T upfc = DjjpfcAco (34) 

D UPF c reduces the generator rotor oscillations. 
Considering (33) and (34), D UPFC can be regarded as 
follows: 

(35) 

1 -K l (s)H(s) 

It can be simplified as follows: 

Dupfc = {K c (s)K 0 (s) + 

Dupfc Kl (s))H(s) 

= F(s)H(s) 

F(s ) in (36) varies with power system operation 
conditions and is defined as follows: 
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F(s) = K c (s) K 0 (s) + 

^UPFC^L is) (37) 

The operation conditions of the power system, 
applying UPFC structure presented in [Sadigh (2010)] 
as a damping controller, depend on the generator 
operation point, power system condition, and the 
UPFC input signal. Considering this fact, the selection 
of UPFC input signals and the most proper of them 
should be in away that the system stability against 
electromechanical oscillations is bonded and it achieves 
its stability as soon as possible. Also, the system 
should maintain its stability as electromechanical 
oscillations occurwhile the damping speed in a power 
system is highly important. As the damping time is 
long, the power system would get damaged. 

UPFC Damping Controller Design 

If it is assumed that a power system has the operation 
domain of v F( / u) , F(s) would be a function of 
operation condition of p and UPFC input signal and 
the controller is designed and installed on power 
system on this base. The desired controller for 
oscillation damping is the one that can guarantee the 
system stability, can acceptably react against the 
generator rotor oscillations, and can damp such 
oscillations in the shortest time. Proper input signal 
selection is one of the most important steps in the 
design of an appropriate controller. For a good 
controller designing, several approaches have been 
evaluated. In this paper, the phase compensation 
method is applied [Chandrakar (2010)]. For proper 
control signal selection of UPFC damping controller, 
three techniques presented in [Wang (1997), Shahir 
(2012)] are applied in the following subsections. 

1) First Technique 

The parameter v F(/v) based proper control signal 
selection for UPFC damping controller design is as 
follows: 

A selected = mil U F i S ’ ^ 11 k ) / ,E ' P G“) ( 38 ) 

Here, the selection aims to making it possible to 
design UPFC damping controller in a way that the 
system would be able to maintain its stability 
against rotor oscillations under any circumstance. 
This is accomplished targeting proper control 
signal selection in order to obtain the power system 
stability in minimum value of V F(//) . 

2) Second Technique 

One of the most important issues considered in the 
designing process of damping controllers is to 
achieve the maximum control domain with 


minimum expense. Therefore, the proper input 
signal is the one that has the maximum effect on 
UPFC damping controller. Proper input signal for 
UPFC damping controller design on the basis of 
achieving minimum control system costs and 
maximum control signal effect is shown in (39): 

F(s, Mselected ,u K ) ^ 
u K £ m E1 , m E 2 , 8 E l , 8 e 2 


3) Third Technique 


In third technique, (40) is applied to select the most 
appropriate input signal considering system 
operation domain, V F(//) . This is done under 
several conditions of x F(/v) to guarantee the power 
system condition and tries to select an input signal 
that would not considerably vary as V F(//) varies 
and has a constant and invariable reaction against 
generator rotor oscillations in V F(//) domain. If the 
influence of damping is increased, other system 
parameters and its operation condition under 
different 'i'(ju) values might be undesirably 
affected. On the other hand, significant decrease in 
the influence of damping controller can decrease 
system stability level and make the power system 
unreliable. In (40), a technique is introduced to 
achieve power system stability against rotor 
oscillations with the least impact on other power 
system parameters and components. However, in 
the cases, the results are not meaningful or not 
expectable; it is possible to apply (39) to guarantee 
the results obtained from (40) for input signal 
selection of damping controller designing. 


^ selected 


max^ F(s,ju, u k ) - min ; , F(s,ju, u k ) 
min fl F(s,jU,u k ) 


(40) 


UPFC Controller Design 


As shown in Fig. 3, the UPFC damping controller 
consists of H(s) feedback transfer function. In [Shahir 
(2012)], the H(s) function is investigated in conventional 
UPFCs damping controllers applied in single-machine 
power systems. In this paper, by developing the 
method in [Shahir (2012)], H(s) is considered as 
follows: 


ov 7) _ sT w K u 1 + sT 2 1 + sT 4 
^ ’ l + sT w l + sT u 1 + sT] l + sT, 
Relation (41) can be expressed as follows: 


(41) 


H(s) = \H\Z0 (42) 

where //| is the amplitude and 0 is the phase angle 
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of H(s ) . Considering (37), F(s) is rewritten as follows: 
F(s) = \F\Z</> (43) 

In (43), |F| is the amplitude and tp is the phase angle 
of F(s) . Both II (s) and F(s) transfer functions 
considerably affect D UPFC and oscillation damping 
capability of power system (see (36)). 

In UPFC damping controller design, the controller 
tries to apply a pure positive torque on investigated 
system using the phase compensation technique. 
Based on (36), F(s) is a function of 'h(//) and it is 
possible to achieve an ideal controller that provides 
positive torque for generator oscillations damping, just 
by controlling H(s) parameters. Considering (36) and 
(42), in order to provide a pure positive torque and to 
design on phase compensation technique, the 
following relation is valid: 

H - ° UPFC AO = A-<j> (44) 

F 

where the amplitude and the phase angle of H(s) is 
obtained and regulated by (44). 

Simulation Results 

One of the most important aims of UPFC installation is 
to control power system against low frequency 
oscillations happening as a distortion occurs in power 
system and the generator rotor is driven out of its 
synchronous mode. UPFC is in charge of proper 
responding as such oscillation occurs. Fig. 4 shows the 
power system illustrated in Fig. 1, where the A a> error 
signal is sampled and applied to UPFC through U 
input signals. In the next parts, the proper signal is 
selected through analytic methods and UPFC 
controller regulation approaches are given. Nonlinear 
simulation is then accomplished on the investigated 
system. 



FIG. 4 UPFC TO SIMB CONNECTION 


F(s) is a function of UPFC damping controller input 
signals (37) and directly relates to 'F ( // ) and varies 
with 'Ft //) and changes AT UPFC magnitude. Therefore, 
F(s) can be applied as a parameter of UPFC controller 
signals evaluation and the most proper control signal 
recognition. The most proper UPFC controller signal is 
selected utilizing techniques introduced in previous 
section and investigating F(s) under different *F( //) 
conditions considered as Table 1 in this paper. The 
magnitudes and the angles shown in Table 2 are 
obtained evaluating F(s) under //, and // 2 operation 
condition applying Nichols [Wang (1997), Shahir 
(2012)] approach. Fig. 5 shows the oscillation 
frequency in V F(//) . It is initially observed that the 
control system is stable for all 'Ff//) values if the control 
system is designed under // 2 operation condition. 
However, this requires further attention especially 
about the UPFC device applied in this paper whose 
operation is controlled by four input signals of A m El , 

Sm E2 , A 8 e x , and A S E 2 ■ In conventional UPFCs applied 
in power systems, the number of input signals is 
decreased to three signals and the investigations are 
carried out on these three inputs due to the incapability 
of using AS b [Wang (1997)], the phase angle deviation 
rate of the series voltage-source voltage. 


TABLE 1 SYSTEM OPERATION CONDITION IN ¥(//) DOMAIN 


operation condition 

V, 

K 

p 

Mi 

1.0/71/ 

1.0/71/ 

0. 1 pu 

Mi 

1.0/71/ 

1.0/71/ 

1.0/71/ 


TABLE 2 THE AMPLITUDE AND PHASE ANGLE OF F(s) UNDER Mi AND 
/J 2 OPERATIN CONDITION 


u k 

FO)| to dB 

phase angle to degree 

A 

Oi 

A 

Mi 

m u 

32.3789 

419.6940“ 

294.8861" 

0.2160 

S Ef 

134.9189 

430.3020“ 

407.4034" 

12.1071 

m E,2 

14.8118 

124.6427" 

313.9723" 

1.4827 

S E,2 

41.7253 

507.4505“ 

300.6187" 

0.3162 


® ® ® 


7.5 


© 


r\(pu) 


Proper Input Signal Selection 


FIG. 5 OSCILLATION FREQUENCY CREATED IN 
INVESTIGATED SYSTEM 
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1) Case study 1 

In this section, the approach mentioned in (38) is 
applied to detect the controller signal. Based on 
(38), the signal that has the minimum value under 
H x operation condition can be selected. Therefore, 
applying (38) and based on the results shown in 
Table 2 results for the most proper signals of 
controller design are m El and S E2 , respectively. 

2) Case study2 

Relation (39) selects the controller signal under ju 2 
operation condition in a way that the input signal 
has the maximum amplitude. Therefore, based on 
the results shown in Table 2, the proper signal for 
controller design is S EA . 

3) Case study 3 

Considering (40) and Table 2, and selecting 
u k = m E l , u k =S El , u k = m E 2 , and u k =S E2 , as 
different input signals, the optimal input signal 
(«, selected ) is obtained using u k = m E ] . 

UPFC Controller Parameters Regulation 

As mentioned before, the UPFC controller parameters 
are regulated applying phase compensation technique 
(relation (44)). The UPFC controller parameters 
displace the Eigen values and are regulated in terms of 
UPFC damping controller input signals in order to 
create desired D UPFC in 'IT//) for generator rotor 
oscillations damping and maintain system stability 
[Shahir (2012)]. As the system Eigen values are less in 
terms of a specific input signal in damping controller, 
the input signal is more effective in damping 
controller and as the UPFC controller gain is less, 
lower K (/ , the UPFC damping controller can be 
designed and utilized with lower cost. The UPFC 
controller parameters applied in this paper are as 
7] = 0.9s , T 3 = 0.9s , T u = 0.01s , and T w = 10s . Other 
UPFC controller parameters are obtained from (44) by 
input signal change. The system Eigen values and 
other UPFC damping controller parameters are 
respectively shown in Table 3 and Table 4 under //, 
and ju 2 in terms of different input signals. Based on 
the results shown in Table 3, the most appropriate 
signals for UPFC controller input under fi x operation 
condition is S E , from effectiveness point of view and 
is m E2 from the view point of cost. Based on the 


results obtained from system analysis in terms of each 
UPFC damping controller signals, the maximum 
damping is achieved selecting S E l as the UPFC 

damping controller input signal while the minimum 
cost is achieved under /u 2 operation condition in 
which m E2 is selecting. 


TABLE 3 UPFC CONTROLLER PARAMETER UNDER /J, 


Input signal 

UPFC controller parameters 

Eigen values 

"Cl 

K u = 0.9, T 2 = 0.2, T 4 =0.15 

s = -0.1302 + y0.0460 

S E , 1 

K u =2.3, T 2 =0.12, r 4 =0.10 

5 = -0.4718 + 70.3120 

m E , 2 

K u = 0.6, r 2 = 0.65, T 4 =0.09 

s = -0.1368 + yO.0101 

S E , 2 

K u = 0.9, r 2 = 0.55, r 4 =0.20 

s = -0.2381 + y'0.1 167 


TABLE 4 UPFC CONTROLLER PARAMETER UNDER jU 2 


Input signal 

UPFC controller parameters 

Eigen values 

"C.i 

K„= 1.1, T 2 =0.95, r 4 =0.95 

,s = -0.0470 + 70.0747 

S B,l 

K u =0.65, r 2 =0.25, P 4 =0.15 

5 = -0.2074 + 70.3264 

m E , 2 

a; = 0.3, r 2 = o.72, r 4 = o.65 

s = —0.0727 i 70.0241 

s E , 2 

a:,, =3.3, 7, =0.65, r 4 =0.45 

^ = -0.0903 ± 70.0003 


Nonlinear Analysis 

In this section, the system shown in Fig. 1 is simulated 
under relations mentioned in (l)-(6). The aim is to 
investigate the way system responds to sudden load 
variation caused by distortions. The simulation is 
carried out under the circumstance that a three-phase 
to ground fault occurs in the system shown in Fig. 1. 
The fault occurs in t = 1 sec . In this simulation, the 
operation condition of the system is as Table 1. Other 
simulation parameters are detailed in Appendix C. Fig. 
6(a) shows the simulation results under //, operation 
condition. The results are evaluated for all four input 
signals. The system condition during distortion occurrence 
without applying UPFC damping controller is shown 
in Fig. 6(a). 



(a) 
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(b) 

FIG. 6 SIMULATION RESULTS (a) UNDER u, (b) UNDER /u 2 


Under such condition, the system starts to oscillate 
and the stability is weakened. Also, Fig. 6(a) shows the 
system condition under distortion occurrence and at 
the presence of UPFC damping controller designed for 
m E | , S El , m E 2 and S E 2 input signals, respectively. 
According to the results shown in the figure, the 
system stability is reinforced and the oscillations are 
rapidly damped under //, and applying u k = S E , , 

while the system damping is slow under u k = m E 2 

condition, which leads to the defection of the power 
system equipments. Fig. 6(b) shows the system 
response to the created distortion when it operates 
under ju 2 operation condition. In this figure, it is 
illustrated that the investigated system oscillates if no 
UPFC damping controller is applied. Also, Fig. 6(b), 
shows the assessment results of system damping rate 
improve the application of UPFC input signals and the 
damping controller for m E1 , S E l , m E2 and S E 2 , 
respectively. According to the results shown in Fig. 
6(b), the system is more stable against the oscillations 
with u k — § E | . Also, it is not reliable with u k =S E2 

and the oscillations are damped very slowly. 


Conclusions 

Although several approaches are used to select proper 
input signals, different standards can be considered to 
select such signals for UPFC damping controller 
installed in power systems (such as the system 
investigated in this paper). If the purpose of damping 
controller design is just to control generator rotor 
oscillations and the economical concerns are not under 
consideration, the most appropriate UPFC damping 
controller signal is the one applied to system desirably 
in all '?(//) and damps oscillations rapidly and 
controls the system. Therefore, it is possible to utilize 
m E2 as input signal. If the economical concerns and 


costs reduction interests are considered in power 
systems designing process, the proper controller is the 
one that possesses lower UPFC controller gain for 
specific input signal in all v F(/0 domains. Therefore, 
S E l can be the most appropriate option. The stability 
of the controllers designed in power system is 
significantly important. The stability of the designed 
controller can be explained as such statement that the 
damping controller input signal is effective in all 
'T (ju) domains and operates invariably. Considering 
such a standard, selecting m El can guarantee the 
designed controller's stability in all V F(//) domains. 


APPENDIX A 


The constants used in (l)-(6) is as follows: 


XMe-X E ,) 

x,x E3 


/ Cl? — 


x a3 x te A 


a 3 - . v (1 JX 2 X 1 , 3 ); b { - 

J a bb 


XJju 

(x EA -xj(x a3 +x a2 y 


X t X E , 


b 2 = 


c, = 


x te (x ai -x a2 ) 

x t x EJ 


■b-,=- 


X B 


X 2 (X ai +x b2 ) + 


(X a 2 + X bi+ X h 2 ){X EX -X te ) + j X x 

x,x EJ 


X te (X a2 + X b 3 +X b 2 )-jX l 

x,x E , 


X 


E, 2 




jX B 


X, 
X t 


bv 


+JX 2 {x a3 +x b3 +x a2 ) 


d\ = 


d"2_ — 


(X te -X E ,)(X b 3 + X b 2 )-jX 1 


x t x E , 


-X te (X b 3 +X b2 ) + jXf 


d 3 = 


JX B b 


x,x EA 


~jX 2 (X b3 +X b2 ) + 


x k 


X 


E, 1 


APPENDIX B 

The constants used in (11)-(15) are as follows: 
K t =[( 0 , +a 2 )(V tq + X q I td )]V b cosS 
- [(a t + a 2 )(V t d - X' d I t q )]V b sin 8 


K 2 =(V t j -X' d I Eq ){c 3 +c 2 ), K 6 = (q +c 2 ) + 


V, 


K 2 =\ + (X d -X' d )(c x +c 2 ), K 4 = -(ci\ +a 2 )V b sin<? 


K 5 = («j + a 2 )X' d V b sin S + (a [ + a 2 ) 


by ^ 

v t.d 


K 


X v V b cosS 
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3 a,m F 7 V h 

K-j = — 1 (cos 8 sin S E2 - sin 8 cos S E2 ) 

4Q C 


?> a 2 m E,l V b 


(cos8sin8 EX -s'm8cos8 EX ) 


3c 2 m E \ 


- (cos S E2 )+ - ' (cos S El ) 


K 9 =- —— — (f cos S F 2 + f 2 sin S E 2 ) 


- 77— (/ 3 cos S EX + f 4 sin S E l ) 
^10 = - "77^" (sin (5 + cos 8 ) , K u = 


K p ,E,2={V t , q +X q I t4 ) -d—X F„ c sin^, 2 


- («i + a 2 - x ' d I ttq )f / /, sin 8 


W,, q +X q I, d )cos8 EX 


~ ( K,d ~ x 'd l t , q ) sin 8 e 1 ]w £ jF* 

X p,SE,2~ 0 /r t,q +X q I t,d)\^~^~ m E,2^dc Cos ^E 


(K,d X' d I, ) m E,2^dc S i n ^E,2 


^d 2 m E2 V dc 


(cos cos S E 2 + sin , sin i) £ 2 ) 


b x m E2 d x m El 

K u = 1C (cos 8 E1 + sin 8 F 2 ) + ^ ’ (cos^, +sin^,) 


'id 2 m F \V dc 


(cos <5 £1 2 + sin cJ £ ) + W| 


^,£,1 =7^ (cos ^ +sin 
^,£,2 = 77^ (cos j £i2 + sin S E2 ) 

^se 

d]Tn F ]Vjr 

K s,se,\ = ' (cos 8 EX - sin S EX ) 

^ se 

K s,se,i = h ' m ^ VdC (cos S E 2 - sin S E 2 ) 

Kqd = f^y^ m E,2 cos S E2 J dl+ ^ d A m EX cos 


K pd ={V, d -X' d I, q ) 


~~~\ m E , 2 cos 8 E2 + 


dl \ dl j «£,1 cos <? £>1 + (F, ? + X q I tyd ) 

■ c- r ^1 ~t d 2 ''i ■ c- 

— ^ I w £,2 sm S E 2 + 1 — — I m E X sin S E X 


v (^if^i+^ 2^1 • c jd x +d 2 '\ • 

K vd = sm 8 e , 2 + — - — J m El Sin^ a 


(b x +bA s 

X q~ 2 — w! £,2 COS £>£ >2 + 


d\ + 


m E \ cos S E X X d 


+ ^2 


Kic cos 8 e x , 


K q,E ,2 = j (rfc COS ^ £i2 

X q, 5 E,\ = ~[ ^ i+ 2 2 ^ m E,l V dc Si 11 S E,\ 
X q,SE ,2 ~ _ f * 2 " \ m E, 2 ^dc si n ^£,2 


W t , d -X' d I t , q )^8 Ed 


+ (K,q + X qdf,d ) Sin S E X ]V dc 


3b 2 m E,2 V dc 


(cos 8 e 2 + sin 8 e 2 ) 


3 b 2 m EX V dc 

+ (cos S E 2 cos S EX + sin S E 2 sin S E X ) + w 2 


3d 3 m E2 m EA V dc 


3d 2 m E i V dc 


(- sin S EX cos S E2 + cos S E X sin S F 2 ) 


+ — 1 — ; (-sin<y £1 cos<y £1 +cos<J £1 sin<y £1 ) + w 3 

8 C dc 


3 b 3 m E 2 V dc 


(-sm^ £ 2 cos S E 2 +cos S E 2 sin S E 2 ) 


3b 2 m El m E2 V dc 


(-sin 8 e 2 cos S EX + cos S E 2 sin S E , ) + w 4 




X' d cos 8 E ^V dl 


K v ,e ,2 = [ 1 y~)[ h^AxqVdc sm S Ea - f j X' d 


d^dc C °S 8 E 2 


d x + d 2 


7 “ \ x 'd m E,i^dc sin^ej + — X q m 6 

z J \ "t J 


_ V t,dYb x +b 2 ') v r/ ^ 
v,SE,2 ~ ^ 2 j X q m ^,2Vdc COS ^7 


b\+b 2 \ v , T . . r 

X d m E, 2V dc sm d E 2 


/,- |)»„co s3 „ + |L w co s3ij 


fi= -7- m £,2 sin 8 e2 + -X m EX sin S EX 


fi - "7" m £,2 COS<5 £ 2 + -X OTg jCOSt^j 
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f- |V 

4 - „ 


w, £,2 sin S E> 2 + ( ~~ 1 m E ,\ sin S Et \ 


J 


V4 C dcJ 
f 3 N 




C0S &E, 1 + . „ ^£,1, ? s * n ^£,1 

4L * 


W, = 


V4 C dcJ 
3w £ ,i 

V 4Q C , 
/ 3/ «£,2 A 

V 4Q C , 


^£,2,rf cos ^E,2 ' 


I E,\,q C0S ^£,1 


I E,2,q C0SJ £ , 2 


^ 3 ^ 
v4Q c y 
/ 3 '«£,1 A 
V 4Q C , 


/ 3 '«£,2 A 

V 4Q c , 


^ E,2,q S * n ^£,2 


^£,l,rf S * n ^£,1 


^E.l.d s ' n &E, 2 


APPENDIX C 

The investigated power system parameters used in 

simulation process are as follows: 

S E X = 30°; S E2 = -25°; X E l = X E2 = 0.2 pu ; 

X se = X te = °- l P u ; ^*v = X d = i-O^M ; X q = 0.6 pu ; 

xj = blpu ; Q c = c, e = 3.0/?zr ; zn £jl = m E1 = 1 ; 

F rfc = E' q = 2.0 pu ■ V b = 1.0 pu ; M = 2H = 8.0s ; 

./; = 50f/z ; r rf0 ' = 5.044s ; = 0.05s ; K a = 50 ; D m = 0 

REFERENCES 

Chandrakar V.K., Dhurvey S.N., and Suke S.C., "Performance 
comparison of SVC with POD and PSS for damping of 
power system oscillations," in Proc. ICETET, 2010, pp. 
247-252. 

Fang W. and Ngan H. W., "Enhancing small signal power 
system stability by coordinating unified power flow 
controller with power system stabilizer," Electr. Power 
Syst. Res., Vol. 65, No. 2, 2003, pp. 91-99. 

Higorani N.G. and Gyugyi L., "Understanding FACTS: 
concepts and technology of flexible ac transmission 
systems," IEEE Press, New Jersey, 1999. 

Kundor P., "Power system stability and control," McGraw- 
Hill Press, New York, 1994. 

Larsen E.V., Sanchez-Gasca J.J., and Chow J.H., "Concept for 
design of FACTS controllers to damp power swings," 
IEEE Trans. Power Syst., Vol. 10, No. 2, 1995, pp. 948-956. 

Sadigh A.K., Hagh M.T., and Sabahi M., "Unified power 
flow controller based on two shunt converters and a 
series capacitor", Electr. Power Syst. Res., Vol. 80, No. 12, 
2010, pp. 1511-1519. 

Saravana llango G., Nagamani C., Sai A.V.S.S.R., and 
Aravindan D., "Control algorithms for control of real 


and reactive power flows and power oscillation damping 
using UPFC", Electr. Power Syst. Res., Vol. 79, No. 4, 
2009, pp. 595-605. 

Shahir F.M., Babaei E., Ranjbar S., and Torabzad S., "New 
Control Methods for Matrix Converter based UPFC 
under Unbalanced Load," in Proc. IICPE, 2012. 

Shahir F.M. and Babaei E., "Evaluation of Power System 
Stability by UPFC via Two Shunt Voltage-Source Converters 
and a Series Capacitor," in Proc. of IEEE ICEE, 2012., pp. 
318-323. 

Shahir F.M., Babaei E., Ranjbar S., and Torabzad S., 

"Dynamic Modeling of UPFC based on Indirect Matrix 

Converter," in Proc. IICPE, 2012. 

Shahir F.M. and Babaei E., "Dynamic Modeling of UPFC by 
Two Shunt Voltage-Source Converters and a Series 
Capacitor," in Proc. ICECT, 2012, pp. 554-558. 

Shahir F.M. and Babaei E., "Evaluating the Dynamic 

Stability of Power System Using UPFC based on Indirect 
Matrix Converter," in Proc. ICECT, 2012, pp. 548-553. 

Shahir F.M. and Babaei E., "Evaluating the Dynamic 

Stability of Power System Using UPFC based on Indirect 
Matrix Converter," Journals of the Engineering and 
Technology, 2013. 

Shahir F.M. and Babaei E., "Dynamic Modeling of UPFC by 
Two Shunt Voltage-Source Converters and a Series 
Capacitor," International Journal of Computer and 
Electrical Engineering, 2013. 

Song S.H., Lim J., and Moon S.II., "Installation and operation 
of FACTS devices for enhancing steady-state security," 
Electr. Power Syst. Res., Vol. 70, No. 1, 2004, pp. 7-15. 

Sreenivasachar K., Jayaram S., and Salama M.M.A., 
"Dynamic stability improvement of multi-machine 
power system with UPFC," Electr. Power Syst. Res., Vol. 
55, No. 1, 2000, pp. 27-37. 

Wang H.F. and Swift F.J., "FACTS-based stabilizer designed 
by the phase compensation method. I. Single-machine 
infinite-bus power systems. II. Multi-machine power 
systems," in Proc. APSCOM, 1997, pp. 638-649. 

Xia Jiang Chow J.H., Edris A-A., Fardanesh B., and 
Uzunovic E., "Transfer path stability enhancement by 
voltage-sourced converter-based FACTS controller," 
IEEE Trans. Power Del., Vol. 25, No. 2, 2010, pp. 1019- 
1025. 


114 



Electrical Engineering Research Vol. 1 Iss. 4, October 2013 


Yassami H., Darabi A., and Rafiei S.M.R., "Power system 
stabilizer design using Strength Pareto multi-objective 
optimization approach," Electr. Power Syst. Res. Vol. 80, 
No. 7, 2010, pp. 838-846. 

Farzad mohammadzadeh Shahir was born in Tabriz, Iran, 
in 1985. He received the B.S. degree in electronic engineering 
from Islamic Azad University, Mianeh, Iran, the M.S. degree 
in electrical engineering from Islamic Azad University, Ahar, 
Iran, in 2008 and 2011, respectively. In 2004, he joined the 
Iran Tractor Manufacturing Company, Tabriz, Iran, where 
he has been an electrical engineering. His current research 
interests include power system dynamic and power 
electronic converters. 


Ebrahim Babaei was born in Ahar, Iran in 1970. He received 
his B.S. and M.S. in Electrical Engineering from the 
Department of Engineering, University of Tabriz, Tabriz, 
Iran, in 1992 and 2001, respectively, graduated with first 
class honors. He received his Ph.D. in Electrical Engineering 
from the Department of Electrical and Computer 
Engineering, University of Tabriz, Tabriz, Iran, in 2007. In 
2004, he joined the Faculty of Electrical and Computer 
Engineering, University of Tabriz. He was an Assistant 
Professor from 2007 to 2011 and has been an Associate 
Professor since 2011. He is the author of more than 210 
journal and conference papers. His current research interests 
include the analysis and control of power electronic 
converters, matrix converters, multilevel converters, FACTS 
devices, power system transients, and power system 
dynamics. 


115 



